A general control variate method for multi-dimensional SDEs: An application to multi-asset options under local stochastic volatility with jumps models in finance

This paper presents a new control variate method for general multi-dimensional stochastic differential equations (SDEs) including jumps in order to reduce the variance of Monte Carlo method. Our control variate method is based on an asymptotic expansion technique, and does not require an explicit characteristic function of SDEs. This is an extension of previous researches using asymptotic expansions to obtain the control variates for such general models. Moreover, in our control variate method, the regression estimators can be chosen for each number of jump times with a stratified sampling, and improve the efficiency of the variance reduction. This paper also provides the asymptotic bias and variance of our method in terms of its terminal time and a small noise parameter used in an asymptotic expansion method. For an application of our method, we evaluate multi-asset options whose payoffs are expressed as linear functions of the underlying asset price processes in general local stochastic volatility with jumps models, and show a calculation scheme of control variates for Greeks. In numerical experiments, we apply the new control variate method to pricing basket, spread, and average options and Delta of basket options under a ZABR-type local stochastic volatility model with jumps, and confirm that our method works very well.